The lightcone of Godel-like spacetimes 
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Abstract. A study of the lightcone of the Godel universe is extended to the so- 
called Godel-like spacetimes. This family of highly symmetric 4-D Lorentzian spaces is 
defined by metrics of the form ds 2 — — (dt + H(x)dy) 2 +D 2 (x)dy 2 + dx 2 +dz 2 , together 
with the requirement of spacetime homogeneity, and includes the Godel metric. The 



• quasi-periodic refocussing of cone generators with startling lens properties, discovered 



by Ozsvath and Schiicking for the lightcone of a plane gravitational wave and also found 
£NJ ' in the Godel universe, is a feature of the whole Godel family. We discuss geometrical 

properties of caustics and show that (a) the focal surfaces are two-dimensional null 
' surfaces generated by non-geodesic null curves and (b) intrinsic differential invariants 

of the cone attain finite values at caustic subsets. 

bJQ, 

PACS numbers: 02.40.Xx, 04.20.-q, 02.40.Hw, 04.20.Jb, 04.20.Gz 



^ 1. Introduction 

<N ' 

The study of null congruences, null hypersurfaces and in particular lightcones in general 
relativity is complicated by the existence of caustics, i.e. points of intersection of the 
generating null geodesies. Since light rays become focussed in the presence of matter as 
well as shear, caustics occur inevitably in realistic situations as frequently encountered 
in astrophysics. The corresponding strong lens effect is an important astrophysical tool 
[48J. In numerical relativity caustics are less welcome, they act as barrier for current 
characteristic codes [HI [131 E2]. The differential geometry of caustics in a spacetime 
setting is still not well developed, contrary to their mathematical classification using 
methods of singularity theory [3]. This may be due to the fairly complicated structure 
of these objects, involving crossings and singularities. Important steps have been taken, 
among others, by Friedrich and Stewart [19] and by the Newman school [20j-[22j, [27J. 

A way towards a better understanding is the study of curved spacetimes with 
analytically known focal surfaces. Ozsvath and Schiicking presented already in 1962 an 
exact and detailed analytical picture of the lightcone of a plane gravitational wave [32J . 
They found a cyclic structure of the focal set, produced by a semi-periodic re-focussing 
of light rays. An often reproduced illustration of a similar light cone drawn by Penrose 
[37] served as starting point for investigations in global Lorentzian geometry |18j . [6]. 
A very similar focal structure is present on the lightcone of a quite different spacetime, 
the rotating Godel universe [281 HI EES]. In view of the T-duality of higher-dimensional 
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supersymmetric versions of the Godel metric and pp waves [5J [25] it is perhaps not 
surprising that the same type of caustic is present. 

For further analytical studies of caustics it makes sense to discuss spaces of 
high symmetry first, since here the geodesic equations can be integrated completely. 
Furthermore, if spacetime homogeneity applies, all lightcones have the same intrinsic 
geometry, independent of the vertex location. To this class of spacetimes belong metrics 
of the Ozsvath class III (33], which include the Godel metric. They have been studied by 
Rosquist et al [29] . Other examples are the spacetime- homogeneous Godel-like or Godel 
type metrics [H|,[39],[TT],[l0]-[45j,[51J, [9], also generalizations of the Godel metric and 
admitting at least a Killing symmetry. Their lightcone, and in particular the focal 
subset, is the subject of the present paper. 

Section 2 shortly reviews the two-parameter family of Godel-like metrics. Basic 
geometrical properties of these metrics depend on a dimensionless parameter k, 
measuring the influence of rotation on the spacetime geometry, k 2 may range from 
— oo to oo, but in this article we confine the discussion to a range of positive k 2 . The 
fc-sequence coincides with the family of (2 + 1)— dimensional geometries investigated 
by Rooman and Spindel [10], if one flat space dimension is added to Rooman-Spindel. 
Their parameter /i is our k. 

The lightcone geometry of the Godel family with 1 < k 2 < oo (k 2 = 2 corresponds 
to the Godel cosmos [23], [21], [20], [5], [31], [35]) is studied in section 3, based on a 
paper by Calvao, Soares and Tiomno [10] . Further sections consider briefly some limiting 
cases. In section 4 we treat the special case k 2 = 1, known as Rebougas-Tiomno metric 
|40} |4"4"]. Its subspace z = const is the three-dimensional anti-deSitter space AdS^. The 
causal family < k 2 < 1 without closed timelike curves is omitted here, only the static 
degeneration k 2 — > with vanishing rotation is shortly considered in section 5. The 
concluding section notes that a cyclic behaviour of caustics on many lightcones (and on 
null hypersurfaces in general) may be expected as a consequence of the Sachs equations 
[47] for divergence and shear of the generator congruence. 

2. Godel-like metrics 

Raychaudhuri and Guha Thakurta [39] have introduced as " homogeneous spacetimes of 
the Godel type" the metrics 



together with the additional condition of spacetime homogeneity. Spacetime 
homogeneity requires at least one further Killing field, apart from the three translational 
Killing vectors along the axes, which evidently exist. This leads to the necessary 
conditions 



with two parameter Q,l. We replace Q by the parameter k = 2Q/1 and consider 
the sequence of metrics labeled by k. Only real k are taken here and I is assumed 



ds 2 = -(dt + H(x)dy) 2 + D 2 {x)dy 2 + dx 2 + dz 2 , 



(1) 



H'/D = const = 2tt, D" / D = const = I 2 



(2) 
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non-negative. Rebougas and Tiomno found that the conditions (j2J) are also sufficient 
for spacetime homogeneity [ID]. Moreover, it was shown in [1DJ that for the metrics 
satisfying (CQ) and (J2]) a further Killing vector exists, leading to a G 5 group of motions. 
A re-examination of the symmetries by Teixeira, Rebougas and Aman [51] , who dropped 
the more or less implicit assumption of time-independent Killing fields made so far, has 
shown that in the special case k 2 = 1 the group is Gy, the maximal symmetry group 
within the Godel-like class of spacetimes. 

We write the metric (CQ) in cylindrical coordinates (t, r, 0, z) as 

ds * = , dt + 2k sinh2 (Zr/2)rf0) 2 + sinh*(fr) a + ^ 2 + ^ 2 
Z Z 2 
The numbers {k 2 ,l) in the two-parameter family ([3]) specify a metric uniquely: Members 
with different pairs (k 2 , Z) represent different spacetimes. In the limit k — > oo and Z — > 
such that kl = 2Q remains finite, the function sinh (Zr)/Z can be replaced by r, and 
([3]) becomes the Som-Raychaudhuri metric [10]. For k 2 — > 2 one recovers the Godel 
metric and for k 2 = 1 the already mentioned G7 metric is obtained, studied in detail by 
Rebougas and Tiomno [40J. 

For convenience we note some properties of the metric (J3]). The non- vanishing 
components of the Ricci tensor are (we follow the conventions of [50J) 

RO = -k 2 l 2 /2, R, 1 = R 2 = l 2 (k 2 - 2)/2, R 2 ° = kl(k 2 - 1)(1 - cosh (Zr)). (4) 

The eigenvalues A determined from det\R^ — AS "| = follow as 

Ai = 0, A 2 = -k 2 l 2 /2, A 3 ,4 = l 2 (k 2 - 2)12. (5) 

The Weyltensor, also given in coordinate form, has the non-vanishing components 

C0101 = - C 1313 = (k 2 - l)l 2 /6, Com = -^(fc 2 - 1) cosh (Zr)/6, 

C0202 = (fc 2 - 1) sinh (Zr)/6, C0303 = -(k 2 - l)/ 2 /3, 

C0323 = - kl(k 2 - l)(cosh (Zr) - l)/3, (6) 

C1212 = (k 2 - 1) (cosh (Zr) - l)(A; 2 (cosh (Zr) - 1) + 2 cosh (Zr) + 2)/6, 

C2323 = (k 2 - l)(cosh (Zr) - l)(-2A; 2 (cosh (Zr) - 1) - 1 - cosh (Zr))/6. 

There exist several interpretations of the matter tensor as calculated from the 
Einstein field equations (including a cosmological constant). The k — > 00 limit, the Som- 
Raychaudhuri metric, describes the gravitational field of a homogeneous distribution of 
charged rotating dust. For the Godel family 1 < k 2 < 00, the combination of a perfect 
fluid, a scalar field and a homogeneous source-free electromagnetic field may serve as 
matter [IDJ. A perfect fluid description alone applies only to the Godel metric k 2 = 2, as 
shown by Bampi and Zordan [4J. For more recent discussions of Godel-like spacetimes 
in gravity theories derived from Lagrangians which are arbitrary functions of curvature 
invariants, see [12] . |43j . 

The interpretation of the Godel family as solutions of the Einstein or other field 
equations may be considered as dubious in the sense that unusual or unphysical forms 
of matter are involved. Therefore these metrics (except Godel) are not treated in 
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the standard book on exact solutions [SD]- But this aspect is not important for the 
geometrical discussion in this paper. The metrics are mainly interesting for their high 
degree of symmetry. All admit at least a G§ group of motions. 



3. The Godel family 1 < k 2 < oo 

3.1. Null geodesies 

For generic Godel-like metrics, use is made of the results by Calvao, Soares and Tiomno 
(CST) [10], also largely keeping their notation for comparison. The authors follow 
a previous paper by Novello, Soares and Tiomno dealing with the Godel metric [30J. 
They give a complete discussion of timelike geodesies and treat also null geodesies. 
We consider only the lightlike case. The high symmetry of the metric allows to write 
down a sufficient number of first integrals for the geodesic equations Dx^/Ds 2 = 0, 
using the fact that for a Killing field k^, ^-k^ is constant along a geodesic. With the 
Killing translations d t , d$, d z and the corresponding integration constants p t , p<j>, p z or 
equivalently p t , (3 = p z /pu 7 = P^/Pt, three first integrals may be written 

i/ Pt = 1 + kl/2 - A; 2 sinh 2 (/r/2)/cosh 2 (lr/2), (7) 

4/pt = kl/{2 cosh 2 (lr/2)) - l 2 j/(4 sinh 2 (lr/2) cosh 2 (lr/2)), (8) 

z/pt = -13. (9) 

The dot denotes the derivative with respect to an affine parameter s. A further relation 

ds ds 



follows from ^^f-g^ = 0: 



,2 ,,■! , ,2 / sinh (ir/2) l^y ^ 2 



r / Pt 1/3 r C osh(/r/2) 2sinh(/r/2)cosh(/r/2)J ' 
As shown by CST, it is convenient to introduce instead of r another radial variable £, 
which increases monotonically with r: 

f = sinh 2 (/r/2). (11) 

Equation ffTOj) then becomes 

i 2 ip\ = i 2 v e + i 2 (i ~(3 2 + km - zv/4 (12) 

with 

v = k 2 + (3 2 -l. (13) 

The equations fJ7|)- f|T2l refer to the class of all null geodesies. We are interested in the 
subset forming a single cone, e.g., passing through the origin of the coordinate system, 
t = 0, r = 0, z = 0. This subset is obtained by setting 7 = 0: Expanding the rhs of 
(fT2"|) around r = or £ m l 2 r 2 /4 = 0, one obtains (d^/ds) 2 m — l 4 p 2 ^ 2 /4 < 0, hence no 
geodesies with 7 / can pass the origin. On the other hand, every geodesic with 7 = 
passes the origin. 

With 7 = the first integrals simplify considerably. The geodesic equations can 
be integrated completely, leading to the following parameter representation of the light 
cone with vertex at t — 0, r = 0, z = 0: 
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2k ( k \ 2w{k 2 -\) 
t = — arctan tanw , (14) 

r = — arsinh (— — — \ , (15) 



= 0o + arccos ( V / l^/(^^)) 1 (16) 

* = "W (17) 

with e = 1(— 1) for the future (past) cone. We have introduced a new affine parameter 
w instead of s by 

w = lptVv(s - s )/2 (18) 

(note ^ > 0, since k 2 > 1 is assumed, the case fc 2 = 1 is treated separately), w > 0(< 0) 
corresponds to the future (past) cone. Equation (flEj) differs from the corresponding 
equation (47) - restricted to 7 = - in pU]. Both equations are correct, but refer to 
different initial values (f> . We have replaced the CST equation in order to have = O 
at the origin r = 0. 

The cone generators depend on the two parameter j3 and 0o, which represent a 
possible pair of transversal coordinates for the light rays. It appears more useful to 
introduce (primarily for the past cone, but easily extended to the full cone) the two 
angular coordinates 9, if on the sky of a suitable observer at the vertex. The observer 
is assumed comoving with the cosmic fluid with the timelike velocity vector = 5q (in 
the case of the Godel metric with k 2 = 2) or defined by the normed timelike eigenvector 
of the Ricci tensor in general. It then is not difficult to see (e.g. by using the method 
described in the second Appendix in [16]) that (3 and O are related to the coordinates 
9, ip (polar angle and longitude) on the observer sky by 

/3 = cos#, O = ¥?. (19) 

We note some well-known or easily accessible results. From (|T5|) follows that null 
geodesies from the origin re-converge after reaching (for rays labeled 9) a maximal 
radial extension r$, so we always have 

2 . , / sintf \ , nn . 
r < rg = - arsinh I — = I . (20) 



1 V \/k 2 - sin 2 9 

re is zero for rays along the polar axis and in the opposite (antipode) direction (9 = 0, tt) 
and reaches its largest value r m for equatorial rays (9 = tv/2). The hypersurfaces 
r = const are always timelike, in particular r = r m is the so-called "light cylinder" or 
"optical horizon": Evidently, the spacetime region r > r m cannot be reached by null 
geodesies from the origin. 

At the horizon r = r m the coefficient of dip 2 in ([3]) is zero, thus the 0-coordinate 
lines become closed lightlike (for r > r max , timelike) lines. They are not geodesies, 
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however. A theory of non-geodesic null curves in a Minkowski spacetime was developed 
by Bonnor j7]. His approach translates immediately to curved spacetimes. A calculation 
shows that the closed null curves on the optical horizon are null helices with constant 
Bonnor curvatures k\ = 1, k 2 = 1(1 + k 2 ) / (4A;), k 3 = 0. 



3.2. Focal subsets and inner metric 

The equations (|T4l) -(fT71) supplemented by (|T9l map the intrinsic coordinates (w,9,(p) 
of the lightcone to the spacetime coordinates (t, r, 0, z). The critical points of this map 
are those where the Jacobian matrix has not the maximal rank 3. This happens if close 
cone generators intersect. For the critial or focal points all four sub determinants of the 
Jacobian must vanish simultaneously: 

d(r,<j),z) _ Q d(<j>,z,t) _ Q d(z,t,r) _ Q d{t,r,<j>) _ Q 
d(w,9,ip) ' d(w,8,if) ' d(w,9,ip) ' d(w,9,ip) 

A straightforward calculation of f l2~Tj) leads to the condition f(w, 9) = for the focal set, 
where 

f(w, 9) = k 2 sin w cos 2 9 + (k 2 — i)w cos w sin 2 9. (22) 

Another way to find singularities is to look for higher degeneration of the induced 
lightcone metric. Numbering the inner coordinates as y 1 = w, y 2 = 9, y 3 = <p, the 
intrinsic three-dimensional cone metric is determined by (i, k = 1...3) 

dx^ dx u 

Qyl Qyh 

7ifc is already degenerate of rank 2. A direct calculation shows that the only nonvanishing 
independent components are 

4(/ 2 — 2/A; 2 cos 2 9 sin 3 w + k 2 q cos 2 9 sin 4 w) 

722 = J3-3 ~ 2 , (24) 

rT] 6 cos z w 

4k sin 9 cos 9 sin 2 w(f — q sin w) 

723 = rT^E ' ( 25 ) 

l z r]°' z cos w 



^ = a-,: ^ lk 9w ( 23 ) 



^±u sin v mil uu .„ , s 

733 = • (26) 



4g sin 2 ^ sin 2 w 

7v 

To obtain these compact expressions we have introduced - besides the focal function 
f(w, 9) - a non-negative function q(w, 9): 

q(w, 9) = (k 2 + cos 2 9-1) cos 2 w + k 2 cos 2 9 sin 2 w. (27) 

For later use we note that q is zero on some closed ^-coordinate lines, defined by 9 = 7r/2, 
w = (2n — l)vr/2, n = 0, ±1, ±2, ±3... and arbitrary (p in the range (0, 2tt). 

The determinant of the two-dimensional metric fT24l - fr26l) can be written as square 
of a function h(w, 9): 

722733 - 723 = h2 i ( 28 ) 
h(w, 9) = 4:f(w, 9) sin 9 sin w/(l 2 rf). (29) 
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Higher order degeneration of the metric requires h = and is therefore given by 

(i) the set of focal points f(w, 9) = 0, where neighbouring light rays intersect, 

(ii) the set of points with w = nir, n integer, called "keel" singularities in [T6] . 
where all rays with equal 9 and different if meet in a point on the nth keel, a spacelike 
line of finite length, and 

(iii) the pole rays 9 = 0,tt, resulting from the sin #-f actor, i.e. from the singularity 
of the polar coordinate system. 



3.3. Newman-Penrose coefficients on the cone 

Additionally to the intrinsic metric, the geometry of a null hypersurface may be 
described by some of the Newman-Penrose spin coefficients, mainly by divergence and 
shear and their change along a ray. 

To illustrate this we first shortly consider a fairly known example, a generic null 
hypersurface in a Minkowski spacetime. Here the real divergence and complex shear 
evolve along a given ray according to the Penrose equations [36J 

p = (p + w[a a - p 2 ])/f, 0- = 0o// (30) 
with the focal function 

/ = 1 - 2wp + w 2 (pl - o-qVq). (31) 

Po and o"o depend on the two transversal parameter fixing a ray. From (1301) we have 

P 2 - W\ 2 = (Po - Wa)/f. (32) 
This equation shows that a parabolic point (a point with p 2 = \o~\ 2 ) on a Minkowskian 
ray implies that the whole ray consists of parabolic points - provided the denominator 
/ in fl32|) does not vanish. The denominator vanishes and thus both p and |c| diverge, 
if the affine parameter w takes one of the two values 

w f = l/(p ±\a \). (33) 

At each Wf & focal surface is passed, and the sign of p 2 — \a\ 2 changes. The quotient 
j = p/\o~\ remains finite, more exactly, j — > ±1 at a focal point. One also notes that a 
focal point can be considered as degenerate parabolic point. 

A similar behaviour of the first-order invariant j at caustics holds for the Ozsvath- 
Schiicking plane wave lightcone pQ and was found in [16] for the lightcone of the Godel 
metric. The difference is only that in both cases one meets an unlimited number of focal 
points if one moves along a ray. It is easy to see that this holds for Godel-like metrics in 
general: If one starts from a cone metric ^ab {A, B always run 2, 3) with h = a/ det\^AB\ 
and w as running (not necessarily affine) parameter on the generating rays, divergence 
and shear amount can be calculated from 

' = -2^ W = t (34) 
Explicitly we find for the metric f[2~^|) - (r2"6"|) 

Q 
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r] /sinw 

It is seen that both p and |cr| diverge at focal points / = 0. The quantity 1/j 2 = \a\ 2 / p 2 
measures the anisotropic part of distance change to neighbouring null geodesies along a 
given ray. It can be written as 

._ 2 ^ ^fi.V — 2<7) sin w cos 2 w Af 2 k 2 cos 2 6 sin 2 w cos 2 w 
^ (/ cos 2w + q sinw) 2 n(f cos 2w + q sinw) 2 

1/j 2 evidently goes to 1 for / — y 0, the same limit is reached at keel singularities w = nir. 
3.4- Geometry of caustics 

The affine parameter w gives rise to a foliation of the cone, but spacelike surfaces 
w = const have no invariant meaning since w is not uniquely determined. There exist 
however invariantly defined two-surfaces on the cone, e.g. the spacelike "zero divergence" 
surfaces. Here p = 0, and from ([35]) one obtains their equation as 

tan2w _ (£; 2 -l)sin 2 # 
2u7~ ~ k 2 - 1 + (k 2 + l)cos 2 #' ^ ' 

Since the rhs is not negative for the metrics considered here, such surfaces can only occur 
at points where tan2w/(2w) is negative or null, that is in the range (2m — 1)tt/2 < 2w < 
mn, m = ±1, ±2, ±3.... 

Other invariantly defined subsets of the cone are the focal surfaces J 7 (described as 
"points of the second kind" in [32J). Their equation / = can be written as 

tanw k 2 — 1 9 n 

= — r ^tan 2 #, 39 

w k 

thus focal surfaces occur at points with (2n — l)ir/2 < w < nir, n = ±1, ±2, ±3..., where 
tanw/w is negative. Contrary to the zero-divergence surfaces, focal surfaces are two- 
dimensional (finite and, as will be argued, non-geodesic) null surfaces: Solving ( 139]) for 9 
and introducing this function 9f(w) in (fH]) - ([T71) . we obtain a parametric representation 
of J 7 . The intrinsic metric of the focal surface follows from (|24]) - (|26]) as 

ds 2 = f ww dw 2 + 2f wip dwd(p + f w d(p 2 (40) 



with 

T sin 2 w(l + T cos 2 w) 3 
* ww = l 2 (k 2 - l) 4 cos 4 w{l + T)\k 2 (l + T) — 1) ' ( ' 

2Tsin 2 w(l + Tcos 2 w) 2 

V " " l 2 (k 2 - 1) 5 /2(1 +T) 5 /2 cos 2 wv ^2( 1 + T ) _ !)' ( > 

4Tsin 2 w(l + T cos 2 w) . . 

and T = -tanw/w. The range of the coordinates w,(p for the nth focal surface is 
< ip < 2n, (2n — l)vr/2 < w < mi. Since f ww f w — f 2 ^ = 0, we have the metric of 
two-dimensional null surfaces, with metric components depending only on w. 
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Figure 1. The first focal surface T\ of the I — 1, k 2 = 2 (Godel) light cone 
is shown as projection into the 3-space t — const using spacetime coordinates 
x = ?-cos0, y = rsin0, z. The cusps at top and bottom lie on the exceptional rays 
and are the intersection points with the first keel. The projected keel is the part of 
the z-axis between top and bottom.- We add a comment of caution. The 3-space 
t = const with the metric tensor is a curved Riemannian V 3 , while the plot must 
use the Euclidean R 3 . Thus distances and angles cannot be represented correctly. In 
particular, distortions can occur if the map between V 3 and R 3 becomes singular, e.g. 
at points where det |<7,&| = 4(£ + (1 — k 2 )£ 2 )/l 2 is zero. This happens at £ = l/(fc 2 — 1), 
corresponding to the "equatorial line" of J\, the intersection of J-\ with the plane 
z = 0. The plot give the wrong impression that this line has nonzero length, while the 
true y 3 -length is zero, cf. (j44|) . 



We consider the first focal surface n = 1 in more detail. T\ is shown in Fig. 1 as 
projection into the 3-space t = const, using [x — r cos 0, y = r sin 0, z) as spacetime 
coordinates. The surface is smooth except at w — n/2 and w = n. At w = ir/2, f w 
tends to zero, f ww to infinity, their product f ww f w (or ) is finite and equal to the 
constant /(Al 4 k 2 (k 2 - 1)). 

We may try to explain this geometrically. As noted above, the sign of the metric 
component in (JTJ decides whether the closed coordinate lines of (f> are spacelike 
or timelike. Calculated on the light cone, g^ becomes 733, and calculated on a focal 
surface on the cone, 733 becomes / w , the 0-coordinate lines are yj-coordinate lines on 
T\. Contrary to g^, 733 and f w cannot become negative, they reach zero only at 
isolated points or lines. Apart from coordinate singularities, the zeros of 733 are found 
at keels w = nir and as zeros of q. q becomes zero only for the (^-coordinate line at 
6 = 7r/2, w = 7r/2 (or w = (2n — l)n/2,n = 0, ±1,2,3... at other focal surfaces J-" n ). 
This particular line is a closed null line, with zero length also from the viewpoint of the 



Godel-like lightcones 



10 



w . / , r-. — o7r cos w sin if ( sm w cos w — w ) , , 

L(w) = I d<pv%? = mu2 1W „„_„„ „ ; x 2 ' (44) 



cone geometry (it is one of the closed null curves on the optical horizon). The other 
^-coordinate lines on T\ (the "parallels" in Fig. 1) are spacelike for 7r/2 < w < n. 
Their total length 

S7r cos w sin 3 w(sin w cos w — w) 
l 2 (k 2 — l)(wcosw — smw) 

increases for w > n/2 from zero to a maximum and declines to zero for w — > ir, when the 
keel is reached. Here T\ shrinks to the two cusp points, and all components f ww , f wtp , f 9V 
vanish. 

Instead of slicing by <^-lines we can represent T\ by the lines y % = (w,9f(w),(p = 
const), the twisted "meridians" in Fig.l. Their tangent vector is ^ = (1, -^-,0), with 
a norm given by (d9 '/ / 1 dw) 2rs ii2 = fww > 0. Thus these lines are spacelike except at the 
end points w = ir/2, w = ir. 

One can construct a focal surface J 7 in a still different way. A two-dimensional null 
surface always admits a foliation by null lines: The equation JabJ 3 = has solutions 
f A different from zero since det\fAB\ — 0. The tangent lines to these directions can be 
taken as w-coordinate lines of a new (u, v) coordinate system on J 7 . In (u, v) coordinates 
the inner metric of a two-dimensional null surface is represented by the normal form 

ds 2 = F(u,v)dv 2 . (45) 

Explicitly, the transformation from (w,ip) to (u,v) is given for J^i by 

u — w, (46) 

f w dw(w — sin w cos w) , 

v = tp- - (47) 

Jw cosw\/smw — w cos w y< k 2 (sin w — wcosw) + wcosw 

and the metric function F(u,v) depends only on u = w. 

4Tsin 2 w(l + T cos 2 w) 
= l 2 (k 2 - 1)(1 + T) 2 ' ( ' 

The w-coordinate lines on J 7 ! have zero lengths and can therefore be denoted as null 
lines, but they are different from the null geodesic generators of the light cone. From 
the four-dimensional viewpoint they are non-geodesic null curves. It should not be 
too difficult to develop a theory of such curves on null hypersurfaces, analogously to 
Bonnor's theory in [7]. 

A way to visualize focal surfaces on the past cone is to locate them on the observer 
sky. We may think of radiation emitted from different parts of the focal surface. If we 
walk down the cone into the past with an increasing affine parameter \w\, after passing 
the first zero-divergence surface at \w\ = 7r/4, the first focal surface (n = 1) starts at 
\w\ = 7r/2. Radiation from caustic points at this epoch would appear to the observer 
as luminous ring along the celestial equator 9 = ir/2. For larger \w\ the focal surface 
radiation comes in as two luminous parallels, moving from the equator towards the 
poles. The pole (9 = 0) and its antipode (9 = tt) are reached for |iw| = it, marking the 
two singular end points of the focal surface (seen as cusps in Fig.l). The cusps are also 
intersections of the focal surface with the two exceptional rays on the past cone (these 
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rays are in the Godel case related to the rotation direction and its antipode direction 
[16], and present also for k 2 ^ 2). 

Keels (denoted as "points of the first kind" in [32]) are another example of 
invariantly defined subsets on the cone. The keels w = nn, parametrized by 9, are 
pieces of spacelike lines connecting the singular end points of the corresponding focal 
surface. At each keel point labeled with 9 all rays with the same 9 and different (p 
intersect. The light cone metric further degenerates at keels and becomes a matrix 
of rank 1, only 722 = 4/ 2 /(/ 2 ^ 3 ) differs in general from zero. At the two end points, 
where the keel meets the corresponding focal surface, the matrix rank of 7^ is zero, 
all components 7^ vanish here. The keel appears as second vertex in representations 
which suppress the z-coordinate, e.g. in the well-known figure of the Godel cone in the 
Hawking-Ellis monograph [26J. But taken all dimensions into account, the keel is an 
extended spacelike line, which shrinks to a point only in the k 2 — > 1 limit of the Godel 
family. The invariant length of the nth keel is given by 



with E and K as complete elliptic integrals. 

As known from the Godel universe or the lightcone of the Ozsvath-Schucking anti- 
Mach metric [32], focal surfaces, keels and zero- divergence surfaces occur quasi-periodic, 
due to the fact that the focal function f(w,9) is not strictly periodic in w, while all 
other functions are circular functions of the affine parameter. 

3.5. Differential invariants 

Another quantitative description of null hypersurfaces is provided by their intrinsic 
differential invariants. The quantity j defined as quotient of divergence and shear is 
already an invariant, it is the only invariant depending exclusively on the first derivatives 
of the cone metric. A comment is necessary here. While \a\ is always not negative by 
definition, p changes the sign when the ray passes a caustic (focal surface or keel), and 
j is +1 or —1 before and behind this point. Thus our formal definition produces jumps 
in j at these points, as written down without further explanation in [16] for the Godel 
cone. This suggests to redefine the first order invariant as j = Xp/\a\, X 2 = 1, with 
appropriately chosen A = f/\f \ or A = sgn(p 2 — \cr\ 2 ), to ensure that j is a continuous 
function through caustics. For example, for the Minkowski space null hypersurfaces we 
have j — jo+w(l —jo)\cro\ as smooth function of w, while j = p/\cr\ shows the unnatural 
discontinuity at the two focal surfaces. Nevertheless we keep j as abbreviation for p/\cr\. 

Besides j there exist higher-order invariants [33], [15], [31]. For their calculation we 
use a triad formalism [TB] . The degenerate inner metric of the cone can be represented 



where ti is a complex covariant vector intrinsic to the cone. The generator direction e* 




(49) 



by 



7*fc titk -\- titk 



(50) 
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satisfies 7jfce fc =0. To obtain a complete co- and contravariant triad on the cone we add 
further vectors t,t, 7, such that 

Ut = 0, = 1, ^t = 0, 7^ = 1. (51) 

We use adapted inner cone coordinates y 1 = w, y 2 = 9, y 3 = ip with e l = S\. The 
degenerate metric 7^ then reduces to the two-dimensional metric 7ab- Comparison 
with ( I24ll26l) gives (together with t\ = 0, and up to a rotation) 

• / kcosOsinw . 

+ r^(/-? smM )' ( 52 ) 



lr)y/q/2 Irfl 2 cosw^q/2 
t 3 = sin 6>/2q sinw/^r]), (53) 

The contravariant components t are calculated from U = jn~t k the result is (besides 
t 1 = 0) 

t 2 = il V y^/(Af), (54) 



I11 . kl cos 9w2rj , . „. . . 

t = = + % .,. n V '— (qemw-f). 55 

4sin#sinw^g/2 4/ sm^cosw-yg 

Rotation coefficients related to this triad and of relevance here can now be obtained 
from 

p + iv= -t%,i-t% )U (56) 
a = -Pt 2tl -n 3jl , (57) 
T =(t 2 t 3 -t 2 ?)(i 2 ,3-i 3 ,2)- (58) 

One may verify that the expressions for p and |<r| obtained from ( 15611571) agree with (J35]) 
and ( )36|) . As noted, the components of t^t are not uniquely determined. This affects 
some rotation coefficients, but not p,\cr\ and also not the invariants. The freedom could 
(but will not here) be used to reach, e.g., v = in fl56|) . For our choice of the triad the 
real and imaginary part of the complex shear is given by 

2r ] cos 2 w- q (l + 2cos 2 w) + g 

q sin 2w 2f cos w 

_ . , k(k 2 — 1) cos 6* sin 2 6 sin 2 w . . 

3m(a) = — - . 60 

A null hypersurface has in general four second-order differential invariants of the 
inner geometry, written as complex quantities I and J and conveniently expresssed in 
terms of rotation coefficients [T4] . The quantity I is linear in the second derivatives of 
the metric, with derivatives only along the generators and, like j, dimensionless: 

' = rV(— )+2^. (61) 

|cr| V p a / \a\ 



Explicitly we find for the real part / 



3 2k{k 2 — 1) cos 9 sin 2 9(r] sinw — f) 
\ a \ J i = J^T2 • ( 62 ) 
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The imaginary part I2 has a more complicated structure: 

\ a \ J 2 = n ,9 • 2 77 o — ; ; — \ ■ ( 63 ) 



A tedious but straightforward calculation shows that 

^0 = VQ 2 sin 2 w(2g — 77), (64) 

i\ = 2 sin w (~ <? 3 + ^Q ,2 (sin 2 w — 3) + r] 2 q(5 — 4 sin 2 w) — 2rf cos 2 wj , (65) 
«2 = 4g 2 — 2qr] cos 2w — n 2 , (66) 
z 3 = -2k 2 cos 2 6> sin w. (67) 

^2 = Dj/p describes the change of the first-order quantity j along the rays. I\ is a 
measure for the rotation of the shear directions (i.e. directions where the distance change 
to neighbouring rays is a maximum or minimum) relative to the generator congruence. 
If Ii is zero (as for null hypersurfaces in a Minkowski or conformally related spacetime), 
the shear directions always point to the same neighbouring null rays if one follows a ray. 

The complex invariant J has the dimension (length) -1 , is nonlinear in the second- 
order derivatives of the inner metric and involves additionally transversal derivatives 
1 14] . J describes changes of the nullsurface geometry in transversal directions, but is 
considerably more complicated than / and will be discussed elsewhere. 

The behaviour of invariants at and in the neighbourhood of focal singularities is 
of interest. While the rotation coefficients p and a show singularities, the invariants 
tend to have finite values. We have already noted j — > ±1 at focal points and keels. 
Expanding / near / = in powers of / leads to 

16k(k 2 — 1) cos 9 sin 2 9 sin w cos 3 w 

q 3 y/rj 
4 cos 2 w(r) — 2q) 



±UfV l rv ±1 U Dill U Dill UU UUB UU „n , „q» , . 

h = — — / + °(A (68) 



q 2 sin w 



-f + o(f 2 ). (69) 



Remarkably, the second-order differential invariants Ii,l2 vanish at focal points. This 
also holds at keels w = nn: Writing w — nn = x, one obtains for small x 

h = 16(-l) n+1 k{k 2 - l)rr 3/2 cos 9 sin 2 9 x 3 + o{x A ), (70) 
4(-l) n r7 x , 9x 

mr(k z — 1) sin 9 

For comparison we note that null hypersurfaces in a Minkowski spacetime satisfy 
J 2 = 1/j — j, thus J 2 vanishes at caustics, I\ is already zero everywhere. 



3.6. Comments on the Godel case as treated in fT6f 

For the Godel cone (k 2 = 2) some differential invariants have been calculated already in 
[T6] . The present paper uses different four- dimensional coordinates as well as different 
transversal light cone coordinates y A . The latter is motivated by the topological fact that 
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no coordinate system can cover the whole sphere without singularity. The coordinates 
u, v in [16] avoid a singularity in the direction of the rotation axes, they become singular 
in equator directions instead. The polar angles 9, here avoid the equator singularities 
but show the usual pole singularities. The relation between both systems of transversal 
coordinates is given by cos0 = (1 — u 2 )/(l + u 2 ), sin9 = \/2(v 2 — l)/(v 2 + 1). For 
k 2 = 2, eqn. f )39|) thus becomes the focal equation — tan w/w = {y 2 — l) 2 / (6t> 2 — 1 — t> 4 ), 
eqn. (81) in [16] .- We note a misprint in eqn. (48) of [16J: the denominator should read 
/ 2 +4(l + / x ) instead of f 2 . 

4. The Rebougas-Tiomno G 7 metric k 2 = 1 

The case k 2 = 1 was excluded so far, we treat it separately. Rebougas and Tiomno 
introduced this special case as "the first exact Godel-type solution of Einstein's equations 
describing a completely causal spacetime- homogeneous rotating universe" [10]. The 
lightcone becomes very simple in this model. Since the Weyl tensor vanishes for k 2 = 1 
(see ([6])), the spacetime metric is conformal to the Minkowski spacetime, thus also the 
lightcone metric is conformal to the Minkowski cone metric. One obtains in the limit 
k 2 — )• 1, / — > sin w cos 2 9, r] — > cos 2 9,q — >■ cos 2 9 of preceding formulae: 

4sin 2 u? 4sin 2 w 2 

722 = Iwfl ' 723 = °' 733 = ~p~ tan °- (72) 

The square root h of the determinant |7ab|, 

h = 4sin9 sin 2 w/(l 2 cos 2 9), (73) 

vanishes at the points w = nit (the only other zeros correspond to the coordinate 
singularity). All light rays from the vertex w = meet again at the points w = nn (n 
integer), which are also vertices. Thus every pair of focal surface and keel in the k 2 > 1 
family of metrics has collapsed into a single vertex in the limit k 2 — > 1. The shear of 
the cone vanishes, only the divergence differs from zero: 

p = — cot w. (74) 

p increases from — oo at w = to zero at w = tc/2 and decreases again until — oo at the 
next vertex w = ir. The lightcone belongs to a type of null hypersurfaces characterized 
by P 7^ 0, | cr | = and denoted as "type 5" in the classification of [15]. There exist no 
second-order inner differential invariants for this class. 

The high symmetry of the Rebougas-Tiomno metric is reflected by the existence of 
a symmetry group G7 [51] , see also [S] for further discussions. 

5. Static degeneration k 2 — > 0, / finite 

The limit k — > 0, keeping I finite, requires Q — > 0. It represents the static degeneration 
of the Godel family and has the simple line element 

rfs 2 = _ dt 2 + dr 2 + sinh 2 (Zr) ^ 2 + 

I 1 



Godel-like lightcones 



15 



Teixera, Rebougas and Aman have shown that this metric admits a six-parameter group 
of motions Gq [51] . The only nonvanishing components of the Ricci tensor are R rr = —I 2 
and Rgg = — sinh 2 (Ir) with a constant Ricci scalar R = —I 2 , and the Riemann tensor 
has only one independent nonvanishing component R r e r e = ~ sinh 2 (Zr). Thus the 
metric cannot easily be interpreted as solution of the field equations, it is nevertheless 
interesting geometrically because of its high symmetry. The null geodesies starting at 
the origin r = 0, z = are given by (here w > corresponds to the past cone) 

t — —w, r = w sin 9, 4> = 4>o, z = w cos 9. (76) 

The equations show that the cone generators do not re-converge as generally for the 
Godel family but extend to null infinity as in the Minkowski spacetime (Minkowski is 
included for I — y 0). No focal surface or keel exist for finite w. Nevertheless the cone 
geometry and in particular the asymptotic behaviour of the cone significantly differ from 
Minkowski. The cone metric with the intrinsic coordinates y 1 = w, y 2 = 9,y 3 = ip is 
given by 



sinh (Zu> sin U 

722 = w , 723 = 0, 733 = , (77) 

thus the determinant \jab\ is the square of the function 

w 

h = — sinh (Iw sin 9). (78) 

Hence, apart from the vertex w = and the coordinate singularity on the symmetry 
axis, there exist no further singularities on the cone. Divergence and shear of the rays 
follow as 

1 lsin9 . . . . 

p= — — coth [Iw sin 9), (79) 

— UJ £ 

1 Zsin^ , . „ s , N 

a = a = coth (Zw sin 0). (80) 

The divergence increases from — oo at the vertex w = to —I sin 9/2 for w — > oo, if 
one goes down the past lightcone, thus it always remains negative. The (real) shear 
starts with zero at the vertex, becomes negative for increasing w and reaches the same 
negative limit —I sin 9/2 as the divergence for w — > oo. One also has 

p 2 -\a\ 2 = 1 -^^- coth (Iw sin 9), (81) 
w 

which is always positive, thus the lightcone consists exclusively of elliptic points. 

The second-order invariants I\ and J are zero, but I2 is different from zero and 
given by 

4lwS sin 9 (CS — Iw sin 9) . . 

2 = (IwC sin 9 + S) (IwC sin 9 - S) 2 ( ' 

with S = sinh (Iw sin 9), C = cosh (Iw sin 9). One verifies that asymptotically, for 
w — > 00, j — > — 1, I2 — > 0, which are standard limiting values at caustics. 
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6. Cyclic structure on general null hypersurfaces 

At the end of their pioneering paper [32] , Ozsvath and Schiicking ask for the origin of the 
periodicity structure. A certain answer can be given by going back to the Sachs equations 
[47] , the differential equations governing divergence and shear on a null hypersurface M 
in terms of Ricci and Weyl tensor projections into M: 

Dp = p 2 + aa + u, (83) 

Da = 2(p-iv)a + (84) 

with 

D = p^ = j-,p= -p KU tH v , v = it^p\ a = -Pw&P, r = -t^lH^ (85) 

w-^ppY, V = C^p^VP. (86) 

The null vector p^ is the direction of the cone generators, the complex null vector t M 
spans spacelike directions in M orthogonal to p^ . It should be stressed that uj and i/j 
- in spite of their origin as projections of four-dimensional quantities - depend only on 
the metric 7^, they are objects of the inner geometry of M. The rotation coefficients 
p, z/, a, t defined in (1851) agree with those calculated from (!56l) - (!58l) for the Godel family. 

The Sachs equations are the first equations to be solved on M for a characteristic 
initial value problem based on the Newman-Penrose formalism and starting from N ' . 
Thus in a sense they can be considered as the Einstein field equations in a nutshell, 
being nonlinear and ruling the influence of matter (u) on the nullsurface geometry. The 
Penrose equations (130]) follow as solution of the Sachs equations in the absence of matter 
and for a vanishing Weyl tensor (u = %p = 0). For the lightcone of the Godel family, 
calculation of u and ip gives 

k 2 cos 2 6 s . , 

[k z — sin 

M . IS 2(k 2 - lUm 2 ^ 9 , 

jfam = -!> 1 (- q + 2rj cos 2 w), 88) 

qrj 

_ .,. 4:k(k 2 — 1) cos 9 sin 2 6 cos w sin w 

3m V = — ^ 1 • 89 

The solution of the Sachs equations with these right-hand-sides is given by ( 135]) and 
( 1591I60D . 

From the Sachs equations on can derive a differential equation for an area distance 
r (not to be confused with the radial coordinate r), introduced by Dr = —pr. One 
obtains with Q = |er| 2 + oj the Jacobi equation 

DDr + Qr = 0. (90) 

The caustics are found as zeros of r. For the Godel family from ( 1H7]) u > 0, hence 
Q >0. 
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The existence of cyclic focal features on many null hyper surf aces (not only cones) 
may then be considered as property of the Jacobi equation, or, more concretely, as 
property of the function Q. The linear second-order differential equation ( 190]) belongs to 
the most widely studied equations in applied mathematics. Starting with the classical 
papers by Sturm, Liouville and Kneser in the nineteenth century, there exist many 
theorems which indeed prove a cyclic or oscillatory behaviour (with arbitrarily large 
numbers of zeros of r) for certain functions Q > 0. For Q < the solutions are 
non-oscillatory, but this holds also for some Q > 0, e.g. for the Minkowski space 
caustics. The precise dependence of the oscillation feature on properties of Q is an open 
mathematical problem, see |53j . 
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